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I.  Summary  of  Work  Accomplished 


A.  Introduction 

The  basic  goal  of  the  project  was  to  investigate  the  behavior 
of  Josephson  junction  arrays  using  theoretical  tools  from  dynamical 
systems  theory.  Emphasis  was  Dlaced  on  the  stability  of  coherent 
operation  of  arrays,  and  their  sensitivity  near  the  onset  of 
instabilities.  Good  progress  has  been  made  in  these  areas;  all  of 
this  work  has  appeared  (or  is  scheduled  to  appear)  in  the  published 
literature:  for  a  listing,  see  references  1-7  in  Section  III.  This 
primary  work  is  summarized  in  Section  IB. 

In  addition,  we  were  able  to  make  progress  on  a  number  of 
other  problems,  which  share  mathematical  features  with  the 
Josephson  junction  system.  This  was  possible  because  the  approach 
of  dynamical  systems  theory  extracts  features  which  are  not 
sensitive  to  the  details  of  the  physical  system;  consequently,  these 
behaviors  am  typical  of  a  broader  class  of  systems.  The  results  on 
these  other  topics  also  have  been  published  in  the  technical 
literature:  for  a  listing,  see  references  8-18  in  Section  III.  This 
additional  work  is  summarized  in  Section  1C. 


B.  Josephson  Junction  Arrays. 


For  practical  applications,  it  is  most  desirable  to  operate 
arrays  in  the  in-phase  state,  so  that  each  oscillator  has  precisely 
the  same  waveform,  and  all  oscillators  are  in  perfect  synchrony. 
Consequently,  the  first  order  of  business  is  to  determine  the 
conditions  under  which  stable  in-phase  operation  is  possible.  The 
stability  regions  for  a  variety  of  circuit  configurations  has  been 
mapped  out,  first  via  numerical  simulations^]  and  then  via 
approximate  analytic  calculations^, 3],  In  fact,  we  discovered  that 

for  certain  circuit  designs,  the  dynamics  possessed  a  fundamental 
symmetry  that  prohibits  stable  in-phase  orbits[4,5].  We  conclude 
immediately  that  such  designs  should  be  avoided  fcr  any  practical 


4 


application.  (For  example,  this  applies  to  small  capacitance  "point 
contact"  junctions  shunted  oy  a  purely  resistive  load.) 

The  next  step  taken  was  to  investigate  the  response  of  the  in- 
phase  state  to  random  fluctuations.  For  small  noise  strengths 
(corresponding  to  very  low  temperatures),  the  effects  show  up  in 
two  fundamentally  different  ways:  (i)  as  broadening  the  sharp 
spectral  lines,  and  (ii)  as  additional  broadband  "skirts”  of  Lorentzian 
character[6J.  As  the  stability  boundary  is  approached,  these  effects 
become  increasingly  strong  -  that  is,  the  quality  ot  the  output  is 
degraded.  How  great  is  the  degradation  depends  only  on  the 
symmetry  properties  of  the  nearing  instability  In  particular,  if  the 
instability  corresponds  to  a  symmetry  breaking,  we  find  a 
remarkable  robustness  possessed  by  the  array:  though  each  element 
suffers  large-scale  fluctuations,  the  tota'  voltage  oscillations 
across  the  array  remain  small!  On  the  other  hand,  if  the  instability 
is  symmetry  preserving,  the  fluctuations  become  increasingly  strong 
as  the  bifurcation  is  approached[7]. 

For  use  as  parametric  amplifiers,  we  also  considered  the  array 
response  to  periodic  perturbations  (representing  the  input  signal). 
Here,  we  fou^d  that  substantial  amplification  is  expected  only  in  the 
pase  of  symmetry  preserving  instabilities.  In  fact,  the 
amplification  of  this  extrinsic  signal  is  expected  to  be  greater  than 
the  degradation  (mentioned  above)  due  to  intrinsic  fluctuations  as  a 
function  of  increasing  array  size.  Specifically,  for  an  array 
consisting  of  N  elements,  the  extrinsic  signal  enjoys  a  power  gain 

proportional  to  N2,  while  the  intrinsic  noise  shows  a  gam  in  power 
proportional  to  N. 

Finally,  we  discovered  the  possibility  of  a  new  kind  of  noise 
sensitivity  that  can  affect  an  array  of  N  junctions,  which  is  absent 
for  the  single  junction  problem.  This  phenomenon,  called  attractor 
crowding[8,9],  is  the  result  of  global  dynamical  considerations  (in 
contrast  with  the  local  considerations  described  above).  In 
particular,  if  the  stable  in-phase  orbit  coexists  with  a  particular 
type  or  or  on  caked  the  splay-phase  orbit  (as  is  often  the  case  in  the 
Josephson  circuit s),  then  noise  becomes  increasingly  effective  ai 
destroying  the  coherence  of  the  array.  Direct  simulations  of  the 
circuit  equations  have  verified  this  picture[9] .  This  problem  can  bo 
avoided  by  operating  the  J^seph^r,.-.  yjnrhcr.  ai.ay  in  a  parameter 
rpr-irog  «,vheie  no  stable  splay-phase  orbits  exist. 
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C.  Additional  Work 


Globally  Coupled  Arrays.  The  primary  spin-offs  of  the 
research  on  Josephson  junction  arrays  rest  on  the  structure  of  these 
equations,  in  particular  the  property  of  global  coupling.  This 
property  allow  certain  general  progress,  which  can  also  be  applied 
to  certain  laser  systems[1 0,1 1],  and  also  to  other  systems[12].  The 
laser  system  we  studied  was  a  diode  pumped  Nd:YAG  laser  with  an 
intracavity  frequency-doubling  crystal.  We  were  able  to  predict  the 
conditions  under  which  stable  steady-state  behavior  was 
expected[10],  and  also  the  appearance  of  splay-state  orbits[1 1  ]. 

Both  sets  of  predictions  were  directly  verified  by  experiments 
carried  out  by  R.  Roy's  group  at  Georgia  Tech.  In  another  project,  we 
looked  at  a  room  temerature  electrical  circuit  consisting  of  an  array 
of  p-n  junctions,  which  has  structural  similarities  with  the 
Josephson  junction  array  circuit.  We  tested  whether  certain 
combinatoric  ideas  might  be  useful  in  describing  the  degree  of 
coherence  exhibited  by  the  array.  In  fact,  the  results  were 
surprisingly  encouraging[1 2],  and  we  conclude  that  such  a  statistical 
approach  might  be  valuable  in  handling  cases  where  the  array 
elements  are  not  nearly  identical. 

Stochastic  Resonance.  This  concerns  the  behavior  of  bistable 
dynamical  systems  which  are  driven  by  a  strong  noise  source  and  a 
weak  periodic  modulation,  and  the  possibility  of  improved  sianal-to- 
noise  ratios  as  the  input  noise  is  increased.  Previously,  in 
collaboration  with  R.  Roy,  we  demonstrated  this  effect  in  a  laser 
experiment.  The  experiment  set  off  a  number  of  theoretical  and 
experimental  efforts,  including  the  publication  of  our  own  detailed 
theory  (which  formed  the  second  part  of  B.  McNamara's  Ph.D. 
dissertation)  [1  3], 

Self-Qraarized  Criticality.  This  refers  to  a  phenomenon 
observed  in  model  spatially  extended  dynamical  systems  which  (like 
the  Joseohson  junction  arrays)  have  large  numbers  of  metastable 
states,  where  the  statistical  steady  state  is  characterized  both 
spatially  and  tc~p»<a;'y  oy  power  law  fluctuations.  Considerable 
interest  has  been  generated,  since  this  offers  a  paradigm  in  which 
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the  balance  between  steady  external  stress  and  internal  transport 
results  in  complex  spatial  and  temporal  behavior,  rather  than  a 
simple  time-independent  equilibrium.  We  established  new  examples 
of  cellular  automata  which  exhibit  SOC[14],  including  the  first  case 
of  a  completely  deterministic  system[15].  We  have  also  p.uvided 
some  of  the  first  analytic  results  on  the  subject,  and  deduced  that 
the  "critical  state"  has  the  structure  of  a  neutrally  stable  attractor 
in  the  appropriate  phase  space. 

Classical  Analog  of  Squeezed  Fluctuations.  Experiments  by  M. 
Bocko  (U.  Rochester)  previously  showed  that  a  classical  analog  of 
the  squeezed  fluctuations  observed  in  quantum  optics  can  occur  in  a 
driven  p-n  junction  near  the  onset  of  simple  bifurcations.  We 
developed  a  theo ry [1 6]  showing  that  such  squeezing  is  a  generic 
property  of  a  broad  class  of  nonlinear  systems  which  includes  the 
one  studied  by  Bocko  The  theory  suggests  a  limit  to  the  maximum 
amount  of  squeezing,  in  contrast  to  (presumably  more  fundamental) 
quantum  optical  descriptions. 
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II.  Personnel 


During  the  period  of  the  contract,  several  people  who 
participated  in  the  research  described  received  support  or  partial 
support.  Bruce  McNamara  is  presently  an  Assistant  Professor  at 
Reed  College  in  Portland,  Oregon.  Kwok  Yeung  Tsang  is  now  a 
scientist  at  the  Naval  Research  Laboratory  in  Washington,  D.C.  Each 
of  these  two  individuals  had  filled  the  postdoctoral  position  funded 
by  the  grant.  The  grant  lias  also  been  used  at  times  to  support  two 
graduate  students:  Li-Shi  Luo  is  scheduled  to  receive  his  Ph.D.  by 
the  end  of  the  calendar  year  1991;  Steven  Nichols  completed  his 
Masters  thesis  in  the  summer  of  1990,  and  is  presently  working 
toward  his  Ph.D.  in  our  group. 
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